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Abstract 

Hamiltonian description of gravitational field contained in a spacetime 
region with boundary S being a null-like hypersurface (a wave front) is 
discussed. Complete generating formula for the Hamiltonian dynamics 
(with no surface integrals neglected) is presented. A quasi-local proof of 
the 1-st law of black holes thermodynamics is obtained as a consequence, 
in case when 5 is a non-expanding horizon. The 0-th law and Penrose 
inequalities are discussed from this point of view. 



1 Introduction 

Evolution of gravitational field within a finite tube with a time-like bound- 
ary was considered in ^ and then reformulated in Here, we extend 
this description to the case of a wave front (a three-dimensional subman- 
ifold whose internal metric is degenerate). Restricting our result to the 
special case of wave fronts, namely to non-expanding horizons, we obtain 
the 1-st law of thermodynamics of black holes as a simple consequence. 

Hamiltonian formulation of any field theory needs always an appropri- 
ate control of boundary data. To illustrate this point, consider as an exam- 
ple the linear theory of an elastic string. Field configuration of a string is 
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described by its displacement function: R x [a, b] 3 {t, x) — > if{t, x) G R , 
fulfilling the wave equation, where velocity "c" is a combination of the 
string's proper density (per unit length) and its elasticity coefficient: 

Passing to appropriate time and length units, we may always put c = f . 
The dynamics of the system may be derived from the Lagrangian density 



L = -^^/Jd^\g^-'(d^^)(d^^) = ii^f - iv'f) , (2) 

where fi,u = 0,1 and (x", x^) = {t, x), g^,^ — diag(— f , +f ), "dot" denotes 
the time derivative and "prime" denotes the space derivative. Entire 
information about the dynamics of the string may be encoded in the form 
of the following Lagrangian generating formula (see [7| or |S] for its correct 
symplectic interpretation) : 

5L{tp,d„(fi) = df, {p'^Sip) = {dfip'^) Sip + {df^(fi) . (3) 

It contains the definition of the canonical momenta: 1) the kinetic mo- 
mentum 

dL 

^'■^^ ^a(^^*^ = '^' 

and 2) the stress density 

± 1 , 
" ^9(9^^=-^^^^-^ ' 

together with the Euler-Lagrange equation, obviously equivalent to (0: 

Integrating infinitesimal generating formula @ over the entire string [a, b\ 
we obtain the finite generating formula: 

sJl^J {-kSip + TrSip) +^Tr^ Sip^'' . (4) 

Hamiltonian description of the same dynamics is obtained via Legendre 
transformation between vr and ip, putting TrSp — S{np) — ipS-K: 



-SH 

with 



J {-kSip — ipSn) + I^TT^fJi/jj , (5) 



H := 1^ {Tvp -L) = \j^ {-K^ + iv'f) ■ (6) 

This formal expression acquires a precise, infinitely-dimensional, Hamil- 
tonian meaning: 

5H . sn 
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as soon as the boundary terms in are killed by imposing e.g. the 
Dirichlet boundary conditions, i.e. by restricting ourselves to an infinitely 
dimensional phase space of initial data {ip, tt), defined on [a, b] and fulfilling 
conditions: ^p{a) = A, ip{b) = B. Within this phase space we have 5ip{a) = 
5ip{a) — and equations Q hold. 

Consider now the subspace of static solutions: ir — = Lp. Due to Q, 
these are points where the derivative of the functional H vanishes and the 
Hamiltonian formula Q reduces to the formula for virtual work 

5H = ~ [tt^^v^] , (8) 

But, due to J^, 7i is manifestly convex. This implies that every static 
solution gives the minimal value of the Hamiltonian in the corresponding 
phase space. Due to equation Q and to boundary conditions, such a 
solution is given by: tt = and 1^5(2;) = ^4 + (a; — a) ■ Inserting this 
value into (jSJ we obtain the following "Penrose-like inequality" : 

^<H, (9) 

analogous to the gravitational Penrose inequality relating the energy car- 
ried by Cauchy data outside of a horizon S and the energy of a black hole 
corresponding to the same value of appropriate boundary data on S. 

Instead of controlling the string configuration at the boundary, we may 
control e.g. its stress by applying an appropriate force F. This leads to 
the Neumann control mode 7r"'"(a) = Ficft, T^'^'ib) ~ -Fright, which is again 
a legitime Hamiltonian system 

-Sii = y {-k Sip - if Stt) - ^ip (Jyr-^j '' , (10) 

where H, obtained via Legendre transformation between (p and tt^ at the 
boundary, plays role of a free energy: 

H := H + ' = H - = ^ /'(^' - (^')' - 2^V^") . (11) 

Again, the boundary term in (IIUII vanishes due to Neumann conditions 
and the field dynamics reduces to JTJ . However, the new Hamiltonian lllll 
is obviously non-convex. There is, therefore, no "Penrose-like" inequality 
in this mode: static solutions corresponding to stationary points of the 
Hamiltonian TC are not minimal points of such a "free energy" . 

In p3 dynamics of the gravitational field within a time-like world 
tube S was analyzed in a similar way. For this purpose the so called 
"affine variational principle" was used, where the Lagrangian function 
depends on the Ricci tensor of a spacetime connection F. In this picture, 
the metric tensor g arises only in the Hamiltonian formulation as the 
momentum canonically conjugate to F. Later, it was proved in that the 
Hamiltonian dynamics obtained this way is universal and does not depend 
upon a specific variational formulation we start with (actually, it can be 
derived from field equations only, without any use of variational principles. 
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the existence of them being a consequence of the "reciprocity" of Einstein 
equations - see [3 and [HI)- On the contrary, the Hamiltonian picture is 
very sensitive to the method of controlhng the boundary data. A hst of 
natural control modes, leading to different "quasilocal Hamiltonians" , is 
given in 0. A conjecture about the "true mass", based on an analysis of 
the linearized theory is also formulated there. 

The aim of the present paper is to give a generalization of the above re- 
sults to the case when the boundary 5 is a three-dimensional submanifold 
of spacetime M whose internal three-metric Qat is degenerate. 

2 Dynamics of the gravitational field in- 
side a null hyper surface 

Consider gravitational field dynamics inside a null hypersurface S: 




Parameter s = ±1 labels two possible situations: an expanding or a 
shrinking wave front (if S is a horizon, these correspond to a black hole or a 
white hole case). To simplify notation we use coordinates x'^, = 0, 1, 2, 3, 
adapted to the above situation: = t is constant on a chosen family of 
Cauchy surfaces whereas is constant on the boundary S (this does not 
mean that x^ is null-like everywhere, but only on S). Coordinates x^, 
A = 1,2, are "angular" coordinates on the 2-surface dV = V f] S whose 
topology is assumed to be that of a 2-sphere. Finally, x'^, k — 1,2,3, 
are spatial coordinates on the Cauchy surfaces {x'^ — const.} and x"", 
a = 0,1, 2, are coordinates on the boundary S. 

In paper pi we derive the following formula: 

-Sn=-^ f (P^'Sgki - mSP''') + ^ / {Ma - a5\) 

^ikjg {^^''^^9AB-2(^woSA" ^A^Swa)) , (12) 

where 

H = ^ / e°o + f / = f / a;, (13) 

and P*"' denotes external curvature of the Cauchy surface, written in 
ADM form. Moreover, A — V det qab is a two-dimensional volume form 
and a = -^log\g°°\. The remaining objects are constructed from a null 
field K tangent to S. It is not unique, since fK is also a null field for 
any function / on S. For purposes of the Hamiltonian formula 11211 we 
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always use the normalization compatible with the (3+l)-decomposition 
used here: < K,dx° > = 1. Hence, K = do — W^Oa- The vector-density 
A" = XK"^ = (A, —An'*) is uniquely defined on S. Now we define^ 

Lb ■■= -g{db,\7aK) = -^£Kgab , (14) 

Wa:=-< \7aK,dx' > , (15) 

where gab is the induced (degenerate) metric on S. Denoting by g^^ the 
inverse two-metric, we define the null mean curvature: / — g^^lAB (often 
denoted by - see |S]). 

The volume term in 1131 vanishes due to constraint equations Q'^^ = 
0^ g°o is often denoted by NH + N''Hk (see e.g. O), where H is the 
scalar ( "Hamiltonian" ) constraint and Hk are the vector ("momentum") 
constraints, A'^ and A'^'^ are the lapse and the shift functions. Constraint 
equations _ff = i Hk = imply vanishing of Q^'q. 

In 9 we give two independent proofs of the formula dJ. The first 
one is analogous to the transition from formula @ to formula 0. For 
this purpose we use Einstein equations written analogously to ^ (cf. ^): 

SL^d^in^-'SA;,) , (16) 

where := Y^\/R and ~ T^, - S^^F^^,. 

Integrating 111611 over a volume V and using metric constraints for the 
connection F, we directly prove 112|l . However, an indirect proof is also 
provided, based on a limiting procedure, when a family of time-like 
surfaces tends to a light-like surface 5*. It is shown that the non-degenerate 
formula derived in and gives l|12|l as a limiting case for e — > 0. 

The last term in 11121 may be written in the following way 

—A'^Swa ~ Xw^SwA ~ u^SWa — u^waSX , 

where Wa := Awa and := g^^goB- Denoting k :— n^WA — wo = 
—K°'Wa we finally obtain the following generating formula: 

-Sn = ^ J^^ [P^'Sgk, - guiSP"') + ^ j^^{,X5a ~ aSX) (17) 
+ {^Xl^'^'SgAB +2(^KSX^n''SWA)) . (18) 

Quantity k fulfilling: K'^VaK = kK, is traditionally called a "surface 
gravity" on S. Its value is not an intrinsic property of the surface itself, 
but depends upon a choice of the null field K on S (i.e. the (34-1)- 
decomposition of spacetime). In black hole thermodynamics there is a 
privileged time, compatible with the Killing field of stationary solution 
and normalized to unity at infinity. In this case the above formula pro- 
vides, as will be seen, the so called first law of black hole thermodynamics. 

^This is the (2-1-1 )-decomposition of the extrinsic curvature Q"'biK) defined in 
■^In the presence of matter the volume term equals 5*0 — §'""^'0 a-nd also vanishes due to 
constraint equations. 
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We stress that the symplectic structure of Cauchy data, given by two 
integrals in l|17|l . is invariant with respect to spacetime diffeomorphisms 
(see [SI)- Neglecting the last, surface integral and defining symplectic form 
only by the volume integral destroys this gauge invariance. 



3 Dynamics of gravitational field outside 
the null surface 

Consider now dynamics of the gravitational field outside a wave front S~ . 
We first add an external, timelike (non-degenerate) boundary S"*" and the 
situation is illustrated by the following figure: 



3<0 





where dV^ = V Ci , and dV~ — V f] S~ . Because dV~ enters with 
negative orientation, we have: J^^ = fgy+ — [gy- ■ Integrating again 
Einstein equations written in the form llfcil over V, using techniques de- 
rived in and [S| to handle surface integrals over dV^ and formula 11211 
to handle the surface integrals over dV~ , we obtain: 



StT Jgv+ ^ / StT Jgv- iDTT J Qy + 

+ t|-/ (\l''^SgAB~2(woSA'' -A^Swa)) , (19) 



167r 



av- 



where a is the "hyperbolic angle" between V and , whereas Q"*" is 
the external curvature of S'^ written in the ADM form (cf. |S]). The 
contribution to the total Hamiltonian from the external boundary is 
written here in the form of a "free energy" proposed in 0: 



1 



Sir 



= Q%-i?o, (20) 



where the additive gauge Eo is chosen in such a way that the entire quan- 
tity vanishes if dV^ is a round sphere in a fiat space. The internal contri- 
bution to the energy is given by formula 1131 with dV replaced by dV~ . 
Shifting the external boundary to space infinity: dV'^ oo, the external 
energy T-C^ gives the ADM mass, which we denote by M, whereas the 
remaining surface integrals over dV^ vanish. This way we obtain the fol- 
lowing generating formula for the field dynamics outside of an arbitrary 
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wave front S in an asymptotically flat spacetime: 

+ _[\l'^''5gAB+2{>^5\-n^mA)) ■ (21) 



4 Black hole thermodynamics 

In this Section we apply the above result to the situation, when the wave 
front S~ is a non-expanding horizon, i.e. I = (see |5])- In this case 
the "internal energy" 7i~ vanishes. Moreover, Einstein equations imply 
lAB _ Q |-ggg and the definition of Wa reduces to: WaK = —WaK. We 
obtain the following generating formula for the black hole dynamics 

-5M = I [P'^Sgki - gkiSP"') + ^ / (XSa - dS\) 

IbTT Jy \ / STT Jgy- 

+ J_ /■ (i^sX - u^SWa) , (22) 

StI- Jr3V- ^ ' 

where s = 1 for a white hole, and s = — 1 for a black hole. 

The so called "black hole thermodynamics" consists in analysing pos- 
sible stationary situations. By stationarity we understand the existence 
of a timelike symmetry (Killing) vector field. If such a field exists, we may 
always choose a coordinate system such that the Killing field becomes 
and all the time derivatives (dots) vanish. Hence, formula 12211 reduces 
to: 

SM = -^ f (KSX-n^SWA) ■ (23) 
Stt Jav- ^ ' 

We assumed here that is tangent to S. If this is not the case, we would 
have a one-parameter family of horizons. Such phenomenon corresponds 
to the Kundt's class of metrics (see e.g. The known metrics of this 

class are not asymptotically flat. We do not know whether or not this 
is a universal property and we exclude such a pathology by the above 
assumption. 

We have shown in [5| that there is a canonical affine fibration tt : 5 ^ 
B over a base manifold B, whose topology is assumed to be that of a 
sphere . The affine structure of the fibers is implied by the fact that 
they are null-geodesic lines in M. Identity —2lab ~ £Kgab ~ implies 
that the metrics <; on 5* may be projected onto the base manifold B, which 
acquires a Riemannian two-metric tensor Hab ■ The degenerate metric gab 
on a manifold 5* is simply the pull back of Hab from B to S: g — ■K*h . 

The quantity Wa is not an intrinsic property of the surface itself, but 
depends upon a choice of the null field K on S. Indeed, if K = exTp{—-y)K 
then Wa = Wa + daf- In particular, there are on 5 vector fields K such 
that K'^VaK = and, consequently, k = 0. These are null geodesic fields 
tangent to fibers of tt : S — > B. 

In case of a black hole, there is a privileged field K, compatible with the 
time-like symmetry of the solution, which is normalized to unity at infinity. 
This way the quantities k and wa in formula 12311 become uniquely defined. 
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We have, therefore, two symmetry fields of the metric Qab on S: do 
and K. Due to normahzation chosen above, we have < do — K, dx" >— 0. 
Hence, the field n :— do — K — u^Oa is purely space-like and projects on 
B. Moreover, it is a symmetry field of the Riemannian two- metric Hab- 

Because the conformal structure of Hab is always isomorphic to the 
conformal structure of the unit sphere S^, we are free to choose a co- 
ordinate system in which Hab = fhAB (and Hab denotes the standard 
unit 2-sphere metrics). The field n is, therefore, the symmetry field of 
this conformal structure. Consequently, n belongs to the six-dimensional 
space of conformal fields on the 2-sphere. Using remaining gauge freedom, 
we may choose angular coordinates (x^) = {6,(f>) in such a way that n 
becomes a rotation field on the 2-sphere. This means (cf. or ini) that 
there exists a coordinate system in which the following holds: 

n = -^t eumV d'"^ . (24) 

Here, f2'° are components of a three-dimensional vector called angular 
velocity of the black hole, and are functions on created by re- 
stricting Cartesian coordinates on 'M? to a unit 2-sphere. We can also 
set 2-coordinate axis parallelly to angular velocity vector field. After a 
suitable rotation we have: (f2*) = (0, 0, f^), z — = cos 9, and: 

n = -n^. (25) 

dip 

Inserting this into 1)23^ we obtain 



1 

'8^ 



u^SWa = ^SJ , (26) 



where 



J=J^---^ I . (27) 

I dV- 



Stt 



is the ^-component of the black hole angular momentum. 

Up to now we have used only the symmetry of conformal structure car- 
ried by hAB ■ The symmetry of the metric itself implies that the conformal 
factor / is constant along the field n. This follows from the observation 
that the trace of the Killing equation implies vanishing of divergence of 
the field n: 

= dAWAethcD n'^)=n'^^AethcD dA.f , (28) 

where the fact that n is the symmetry field of the metric h has been used. 
Formula I25I I implies that 9^/ — and the conformal factor / must be a 
function of the variable 9 only. 

It turns out that also its canonical conjugate k may be gauged in such 
a way that it is constant along the field n (see for a proof ).'^ 

This result was obtained locally, or rather ^ttasi-locally - i.e. from the 
analysis of the field on the horizon itself. However, the global theorems 



^In case £7 = 0, quantities ft and / are arbitrary functions on S 
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on the existence of stationary solutions possessing a horizon, imply the 
so called 0-th law of thermodynamics of black holes (see |12|V according 
to which the surface gravity k must be constant along the horizon. But 
A := ^ is the area of the horizon 5*. Taking this into account and using 
()2fi|l , we derive from II23|I the "first law of black holes thermodynamics" : 

-sSM = —kSA + nSJ . (29) 

Contrary to the theory proposed by Wald and Iyer in |15|. the 1-st law 
(I29I I is, in our approach, a simple consequence of the complete Hamilto- 
nian formula 12211 . restricted to the stationary case. As illustrated by an 
example of the string dynamics, where formula ||HJ for virtual work was a 
consequence of the Hamiltonian formula a similar "thermodynamics 
of boundary data" may be expected in any Hamiltonian field theory (see 
e.g. |S] for the corresponding analysis of the Maxwell electrodynamics). 
Also a "Penrose-like" inequality (analogous to @ in the string theory) 
is satisfied as soon as the Hamiltonian is convex. We very much hope 
that the gravitational Penrose inequality can be proved along these lines. 
Preliminary results in this direction, based on the analysis of the field 
Hamiltonian in linearized gravity (see |4j), are promising. 
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